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RESUMEN 

El principal objetivo de este artfeulo es presentar la forma de calculo de la 
curva de capacidad sfsmica resistente, que relaciona el cortante basal con el 
desplazamiento lateral maximo, en una estructura con disipadores de energfa por 
fluencia del material, tipo ADAS o TADAS, modelados como un elemento mas del 
portico, empleando la tecnica del Pushover multimodal. 

La matriz de rigidez de los elementos de la estructura, se hallan por el 
metodo de las dovelas, considerando para los disipadores 5 dovelas y para los 
restantes elementos 3 dovelas, se indica con detalle la forma de calculo de esta 
matriz y del modelo constitutive bilineal que se considera. 

Se analizan dos estructuras, una de hormigon de 4 pisos y otra de acero 
de 6 pisos, en cada una de ellas se considera que tiene disipadores ADAS por un 
lado y TADAS por otro, de tal manera que se obtienen curvas de capacidad 
sfsmica para 4 casos. 

Se indican las curvas de capacidad que se hallan en cada modo de 
vibracion y la resultante que se halla al aplicar el criterio de la combinacion 
cuadratica completa CQC. La curva resultante se compara con la que se obtiene 
con un Pushover monotonico en que se trabaja con el primer modo de vibracion. 

Palabras claves: Pushover monotonico; Pushover multimodal; disipadores ADAS; 
disipadores TADAS. 

ABSTRACT 

The main objective of this article is to present the calculation method of the 
seismic resistant capacity curve, which relates the base shear with maximum 
lateral displacement in a structure with energy dissipation by yielding of the 
material, type ADAS or TADAS, modeled as an element of the frame, using the 
technique of multimodal Pushover. 
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The stiffness matrix elements of the structure, is found by the method of 
segments, considering for the dissipators five segments and for the remaining 
elements 3 segments, is indicated in detail how this matrix and the bilinear 
constitutive model considered are calculated. 

Two structures are considered, one of concrete with 4 floors and other of 
steel with 6 floors, each one of them are considered to have ADAS and TADAS 
dissipators, so seismic capability curves are obtained for 4 cases. 

Capacity curves that are in each vibration mode are found and the 
resultant which is found applying the criterion of complete quadratic combination 
CQC. The resultant curve is compared to that obtained with a monotonic Pushover 
when working with the first vibration mode. 

Keywords: Monotonic Pushover; Multimodal Pushover; ADAS or TADAS Energy 
Dissipators. 



1 INTRODUCCION 

Hay dos temas que se abordan en el presente articulo, el primero de ellos 
tiene que ver con el uso de disipadores de energia por fluencia del material 
denominados ADAS y TADAS, colocados sobre diagonales de acero en forma de 
"V" invertida; con estos dispositivos de control pasivo, se aumenta la rigidez de la 
estructura y lo principal se le confiere amortiguamiento a la misma. (Aguiar et al. 
201 5,1 ; Tsai et al. 1 993, Wittaker et al. 1 989). 

En la parte superior de la figura 1 , se observa un disipador ADAS, que 
tiene la forma de un reloj del tiempo, con un ancho b r en la parte superior; b 2 en la 
parte mas angosta; tiene una alturah; el espesor de cada placa es t; y esta 
compuesto por n placas, generalmente de acero, aunque ultimamente se estan 
trabajando con placas de cobre, aluminio y zing, que tienen mas ductilidad que el 
acero. Heresi (2012). 

En la parte inferior de la figura 1, se muestra un disipador TADAS que 
tiene la forma de un triangulo con un ancho b en la parte superior; altura h; 
espesor de placa t; ancho de placa en la parte inferior b t \ numero de placas n. 

El comportamiento de histeresis de los disipadores ADAS y TADAS es de tipo 
bilineal. Normalmente se han colocado sobre diagonales de acero en forma de "V" 
invertida; de tal manera que durante un sismo, una de las diagonales trabaja a 
traccion y la otra trabaja a compresion. A traccion tiene mejor comportamiento el 
acero que a compresion; en Aguiar et al. (201 5,2) se ilustra el calculo del diagrama 
de histeresis fuerza-desplazamiento, de tipo bilineal, que tienen las diagonales de 
acero. 
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GEOMETRlA DE LA PLACA 




t=espesor de lo placa 





Disipador TADAS 

Figura 1 Disipadores de energia por fluencia el material: ADAS y TADAS 

La geometria de los disipadores ADAS y TADAS, garantiza que cuando 
ingresan al rango no lineal, toda la seccion se va a plastificar porque su curvatura 
es constante; en el caso de los ADAS se tiene doble curvatura, por su forma; y en 
el caso de los TADAS, simple curvatura. Nuevamente en Aguiar et al. (2015,2) se 
ilustra la forma de calculo del diagrama momento-curvatura de estos dispositivos 
mediante un modelo bilineal. 

En la figura 1, se aprecia que para los dos disipadores existe un 
corrimiento relativo A entre el borde superior y el borde inferior. Este 
desplazamiento es el que da origen a la liberacion de energia durante un sismo y 
por ende a incrementar el amortiguamiento. 

El otro tema que se aborda, en el articulo, tiene que ver con la obtencion 
de la curva de capacidad sismica resistente, que relaciona el cortante basal, V, 
con el desplazamiento lateral maximo, D t , mediante la Tecnica del Pushover. Es 
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Figura 2 llustracion de la obtencion de la curva de capacidad sismica resistente 
aplicando la Tecnica del Pushover. Fuente: Aguiar (2002) 

En la parte superior de la figura 2 se observa el comportamiento no lineal 
de una estructura sometida a cargas ciclicas, y en la parte inferior se observa que 
la curva de capacidad sismica resistente que se halla al aplicar la Tecnica del 
Pushover es muy aproximada a la que se halla en la curva de la envolvente de la 
grafica superior. 

Esta semejanza fue la que motivo a trabajar en el tema del Pushover que 
no es mas que la aplicacion de cargas laterales hasta llevar a la estructura a un 
punto que se considera el colapso de la estructura. 

En los primeros trabajos, la aplicacion de la carga lateral se la hizo en 
funcion del primer modo de vibracion, lo cual es bastante bueno para estructuras 
regulares, razon por la cual fue muy criticado este procedimiento por algunos 
investigadores entre ellos Chopra y Goel (2001 , 2002) quienes propusieron 
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trabajar con todos los modos de vibracion en lugar de solo hacerlo con el primer 
modo y es el tema que se desarrolla en este articulo. 

2 MODELAJE DE ELEMENTOS 



2.1 Elemento disipador 

Existen soluciones analiticas, propuestas por Whitaker et al. (1987), Tena 
(1997) para encontrar la matriz de rigidez del elemento disipador sea este TADAS 
o ADAS pero en este articulo se prefiere calcular la matriz de rigidez, dividiendo al 
elemento en 5 dovelas, las mismas que se indican en la figura 3. 

b 

0.9b ; 
0.7b 

0.5b | 
0.3b I 
1 1 0.1b I 

TT 

Figura 3 Dovelas consideradas en elementos disipador. 

Cada una de las dovelas es de seccion contante, cuya base esta indicada 
en la figura 3 y la altura de la seccion transversal es el espesor t, estrictamente es 
nt, siendo nel numero de placas, y la longitud es h/S. Al ser de seccion 
constante es directo el calculo de la matriz de rigidez de cada dovela en 
coordenadas globales. Mora y Aguiar (201 5). 

La numeracion de los grados de libertad se inicia en el borde inferior del 
disipador, luego se continua en el borde superior del disipador y posteriormente en 
los puntos de contacto de las placas. En total se tienen 18 grados de libertad; se 
halla la matriz de rigidez por ensamblaje directo y se condensa a las coordenadas 
de los extremos las mismas que estan indicadas en la figura 4, estas son las 
coordenadas globales del elemento disipador. 
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Figura 4 Coordenadas globales del elemento disipador. 

Al ser las dimensiones del elemento disipador muy pequenas, no es 
importante considerar el efecto de corte 0. Para el modelo constitutive bilineal se 
trabaja con la dovela superior que tiene un ancho b y una altura de la seccion 
transversal igual a: n t. En esta dovela se tiene que el momento de inercia /, el eje 
neutro c, el momento plastico z, y el momento de fluencia M Y , valen: 

Hntf nt (1)y(2) 

12 2 



z = 



fo (n t) 2 , ( 3 ) y ( 4 ) 

My = Z * fy 1 ' ' l ' 



La variable todavia no definida es el limite de fluencia del material f y . Para 
hallar la curvatura de fluencia se la obtiene con la siguiente ecuacion. 



0 Y = 



My_ (5) 



E*I 



Donde E es el modulo de elasticidad del material. La curvatura ultima 0 V 
se halla en funcion de la ductilidad por curvatura Finalmente el momento ultimo 
M u se determina en funcion de a que es la relacion entre la rigidez post fluencia, 
con respecto a la rigidez elastica. 

0 U = Ma 0y M a = M Y + a E I (0u - 0 Y ) ( 6 ) y ( 7 ) 

En la figura 5 se presenta el diagrama momento-curvatura para el 
elemento disipador, para los dos casos en que va a trabajar y son: la parte inferior 
trabaja a traccion y la superior a compresion y viceversa. 



Por ser la curvatura uniforme, en los disipadores se considera que el 
diagrama momento curvatura encontrado, rige para todas las dovelas. 
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Figura 5 Diagrama momento-curvatura para elemento disipador 



2.2 Elemento diagonal 

La diagonal de acero del contraviento, trabaja axialmente, de tal manera 
que la matriz de rigidez del elemento es de 4 por 4. (Aguiar 2014). Para el 
comportamiento no lineal, se tienen dos posibilidades, la una que la diagonal 
trabaje a compresion, en este caso se trabaja con el modelo constitutivo indicado 
en la figura 6 y la otra posibilidad es que la diagonal trabaje a traccion, con el 
modelo que se encuentra en la figura 7. 




Figura 6 Modelo bilineal para elemento diagonal, trabajando a compresion. 

La diferencia entre estos dos modelos es la forma como se obtiene el 
punto de fluencia, para cuando la diagonal trabaja a compresion, la fuerza del 
punto de fluencia es el menor valor entre F Y y F Cr . 

F Y = A*f y (8) 

Donde A es el area de la seccion transversal de la diagonal; F CL es la fuerza critica 
de pandeo, de acuerdo al ASCE 41 (2011), para perfiles tubulares, se debe 
encontrar primero la relacion entre la longitud del elemento L, y el radio de giro r. 
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Figura 7 Modelo bilineal para elemento diagonal, trabajando a traccion. 



- < 108 ( 9 ) 

r 

Si se cumple la relacion f cr = 12 psi; caso contrario se obtiene f cr con la 
siguiente expresion. 

1.4* 10 s (10) 

fcr = r^— (psi) 



Para pasar de (psi) a (T/m 2 ) se debe multiplicar por 703.069. Finalmente, 
la carga critica de pandeo se halla con la siguiente ecuacion. 

F cr =Af cr (11) 

Para el caso de traccion la fuerza de fluencia es simplemente F Y . Para 
hallar el desplazamiento de fluencia, A Y tanto a compresion como traccion se 
divide la fuerza de fluencia para la rigidez axial EA. El desplazamiento ultimo, Ay 
es igual a la ductilidad por el desplazamiento de fluencia; finalmente la fuerza 
ultima, F u se halla en funcion de a pero esta vez relaciona la rigidez axial post 
fluencia con respecto a la rigidez axial elastica. 

2.3 Elemento viga o columna 



En la figura 1 , se muestra el modelo basico denominado sistema 1 , a partir 
del cual se halla la matriz de rigidez del elemento, en coordenadas locales 
(sistema 2) y globales (sistema 3). Este sistema basico es el que se va a utilizar y 
se lo presenta a continuacion primero para cuando el elemento trabaja en el rango 
elastico y posteriormente cuando ingresa al rango no lineal. Para los dos casos se 
puede o no considerar el efecto de corte 0 . 
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Figura 8 Sistema de coordenadas 1 , en que no se consideran desplazamientos 

como cuerpo ngido. 



2.3.1 Rango elastico 



En la figura 8 se ha denominado A 0 ,I 0 , al area e inercia en el centro de 
luz, que es igual a la del nudo inicial y final, para el caso de elementos de seccion 
constantes. En general para el rango elastico sea A, I el area y la inercia a flexion. 
La matriz de rigidez, para el sistema de coordenadas 1, es la siguiente. Aguiar 
(2012, 2014). 



K = 



4EI (1 + 0) 
L (1 + 40) 

2EI (1-20) 
L (1 + 40) 

0 



2EI (1-20) 
L (1 + 40) 

4EI (1 + 0) 
L (1 + 40) 

0 



0 

0 

EA 

IT 



(12) 



Todas las variables han sido ya definidas, restando indicar la ecuacion con 
la cual se halla el efecto de corte 0 



0 = 



3/? E I 
GAL 2 



(13) 



Donde G es el modulo de corte; /? es el factor que mide el efecto de corte, 
para elementos de seccion constante vale 1 .2. 

La matriz de paso del sistema 1 al sistema 2 (coordenadas locales) es 
T ± _ 2 y la matriz de rigidez en el sistema 2 (figura 9) se denomina k 2 y se halla con 
la siguiente ecuacion matricial. 



fe 2 — t\_ 2 k l r 1 _ 2 



(14) 
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T = 
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Figura 9 Sistema 2. Coordenadas locales 



Para pasar del sistema 2 al sistema 3 (coordenadas globales) se utiliza la 
matriz de rotacion T 2 _ 3 y la matriz de rigidez en coordenadas globales fe 3 se halla 
con la siguiente ecuacion pero antes en la figura 1 0 se indica el sistema 3. 



T - 

± 2-3 — 





&3 — ^2-3 


fe 2 T 2 


-3 






cos a 


sena 


0 


0 


0 


0 


sena 


cos« 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


cos a 


sena 


0 


0 


0 


0 


- sena 


cos a 


0 


0 


0 


0 


0 


0 


1 



(16) 



(17) 



Donde a es el angulo que forma el eje del elemento con la horizontal, para 
vigas vale cero grados y para columnas 90 grados. 
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Figura 10 Coordenadas locales y globales de un elemento. 

2.3.2 Rango inelastico 

En la figura 11 se presenta el modelo numerico de calculo, para un 
elemento viga o columna, que trabaja en el rango no lineal, en el sistema de 
coordenadas 1 . Sea L pa la longitud plastica en el nudo inicial (El momento en el 
nudo inicial es mayor que el momento de fluencia M Y ); sea L pb la longitud 'plastica 
en el nudo final (El momento en el nudo final es mayor que el de fluencia M' Y ) y L 0 
es la longitud sin dano. (Mora y Aguiar, 201 5). 




X 

> 



Figura 11 Elemento en coordenadas de sistema reducido 

El modelo de la figura 11, considera la longitud de dano, a diferencia de 
otros modelos de plasticidad concentrada, como el de Giberson (1968). En el 
modelo de la figura 1 1 se considera que el elemento esta compuesto por tres sub 
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elementos, que tienen rigideces a flexion: EI a ,EI 0 ,EI b y rigidez axial: EA a ,EA 0 ,EA b ; 
el subindice a, es para el sub elemento de longitud L pa ; el subindice 0 para el sub 
elemento de longitud L p0 ; y el subindice b para el sub elemento de longitud L pb . En 
cada sub elemento se halla la matriz de rigidez en coordenadas globales, en la 
forma indicada en el sub apartado 2.3.1, considerando o no el efecto de corte 0. 

Se ensambla la matriz de rigidez, teniendo en cuenta los grados de 
libertad indicados en la figura 1 1 , de tal manera que se halla una matriz de 9/9 y 
se condensa esta matriz a las tres primeras coordenadas. (Mora y Aguiar, 2015). 
De esta manera se halla la matriz de rigidez para el sistema 1 , pero considerando 
el comportamiento no lineal de los elementos. En Aguiar et al. (2015,3) se ilustra la 
forma de proceder a partir de este modelo, para cuatro condiciones: i) Todas las 
secciones se encuentran en el rango elastico; ii) Se forma la rotula plastica en el 
nudo inicial y en el final no; iii) se forma la rotula plastica en el nudo final y en el 
inicial no; iv) se forman rotulas plasticas en el inicial y final. 

La relacion momento-rotacion de los elementos, se obtuvo a partir de los 
criterios y tablas, recomendadas por el ASCE 41 (201 1) tanto para elementos de 
hormigon como para acero. En la figura 12 se presenta a la izquierda el diagrama 
momento curvatura con el cual se halla las rigideces a flexion y a la derecha se 
indica el diagrama momento-rotacion. En Aguiar et al. (2015,2) esta muy detallado 
la forma de encontrar: el momento de rotacion en la fluencia; la longitud de la 
zonas plasticas; el calculo del momento de fluencia con el modelo de Y. Park 
(1985) para hormigon armado; el calculo del momento de fluencia para elementos 
de acero; la forma de pasar del diagrama momento rotacion al diagrama momento 
curvatura. 




Figura 12 Relacion momento curvatura y momento rotacion, con modelo ASCE 41 

Un tema muy importante es la forma como se trabaja en la zona residual 
~DE donde el momento residual M r -c M Y . El factor c recomienda el ASCE 41 . De 
tal manera que en esa zona se ha superado el momento ultimo, aqui la 
deformacion del hormigon es mayor a 0.004, no se ha formado la rotula plastica 
pero la seccion todavia tiene capacidad a flexion. En Mora y Aguiar (2015) esta 
descrita la forma de trabajar en esta zona. 
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3 MATRICES DE RIGIDEZ, MASA Y PROPIEDADES 
DINAMICAS 

Posteriormente, se presentara la curva de capacidad sismica resistente, 
mediante un pushover multimodal, de la estructura de 4 pisos de hormigon armado 
con disipadores de energia, indicada en las figuras 13 y 14. Pero en este apartado 
interesa ilustrar cual es el modelo de analisis, sus grados de liberad, las 
coordenadas principales y las coordenadas secundarias, esto se encuentra en la 
figura 15. 



0.75 T/m 




Figura 13 Estructura de hormigon armado con disipadores de energia. 

Es de interes conocer la matriz de rigidez y la matriz de masas, al inicio del 
pushover, asi como sus periodos de vibracion, modos de vibracion y factor de 
participacion modal. Pero a mas de ello que se presente un programa que 
resuelva lo indicado, utilizando la libreria de programas del sistema de 
computacion CEINCI-LAB. 

La curva de capacidad sismica que se presentara posteriormente, 
considera dos casos: el primero que la estructura tiene disipadores de energia 
TADAS, los mismos que se indican a la izquierda de la figura 14 y el segundo que 
la estructura tiene los disipadores de energia ADAS que se indican a la derecha 
de la figura 14. El programa que se indica mas adelante, en este apartado 
considera que la estructura tiene disipadores TADAS. 

Notese que el disipador de la primera planta, tiene 6 placas de acero y es 
diferente del disipador de la segunda y tercera planta, que tiene 4 placas de acero. 
El espesor de las placas es de 1 cm., cantidad con la que espera la plastificacion 
de los disipadores. 
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Figura 14 Disipadores de energia considerados en estructura de 4 pisos. 

Se considera que el disipador es un elemento mas de la estructura, al 
igual que cada una de las diagonales que conforman la "V" invertida. Para 
entender el programa es necesario conocer la numeracion de los nudos y 
elementos, los mismos que se indican a la izquierda de la figura 15. 
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Figura 15 Numeracion de nudos, elementos y grados de libertad. 



A la derecha de la figura 1 5 se muestran los grados de libertad, primero se 
ha numerado las coordenadas principales de la 1 a la 7. Se destaca que el 
desplazamiento lateral 1, en el nudo formado por la union de las diagonales y el 
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disipador, es diferente del desplazamiento lateral del primer piso que se ha 
denominado 2; precisamente esta diferencia de desplazamientos es la que origina 
la disipacion de energia. 



3.1 Matriz de rigidez 

En la grafica de la izquierda de la figura 1 5, se observa que primero se han 
numerado todos los elementos de hormigon armado, que van del 1 al 16. Luego 
las diagonales de acero que van del 17 al 22 y finalmente los disipadores del 23 al 
25. Es importante numerar por separado para poder encontrar la matriz de rigidez 
de la siguiente forma. 

KT = KH + KDIAG + KDIS ( 18 ) 



Donde KH es la contribucion a la matriz de rigidez de los elementos de 
hormigon armado; KDIAG es la contribucion de las diagonales de acero, de los 
montantes, a la matriz de rigidez y KDIS es la contribucion de los elementos 
disipadores a la matriz de rigidez. 

Las secciones de los elementos de hormigon armado se indican a la 
derecha de la figura 13; las columnas son de 50/50 y las vigas de 40/40; la 
armadura longitudinal y transversal de las columnas es la misma en todos los 
elementos; las vigas tienen un armado diferente en los apoyos y en el centro de 
luz pero en todos los pisos el refuerzo es el mismo. 

Las diagonales de acero, estan formadas por perfiles tubulares cuadrados 
de 100/100/10 mm., como se indica en la parte superior derecha de la figura 15 y 
como ya se indico la geometna de los disipadores se halla en la figura 16. 

Cada una de las contribuciones a la matriz de rigidez, de los diferentes 
materiales, se realiza por ensamblaje directo, trabajando con la matriz de rigidez 
de los elementos en coordenadas globales y con los vectores de colocacion. Para 
el ejemplo la matriz de rigidez KT es de 37 por 37. Aguiar (2014). 

Una vez que se halla la matriz de rigidez total, se condensa a las 
coordenadas principales. Cuando se numera primero las coordenadas principales 
hay solo dos formas para la condensacion estatica, la una que involucra el calculo 
de una matriz inversa y la otra la solucion de un sistema de ecuaciones. 

En cambio cuando las coordenadas principales se numeran al ultimo, hay 
tres formas de calcular la matriz de rigidez condensada, las dos primeras han sido 
ya indicadas en el parrafo anterior y la tercera, que es la optima, realizando la 
triangularizacion de la matriz KT hasta la fila n - na. Siendo n el numero total de 
grados de libertad y na el numero de coordenadas principales (Aguiar 2012, 2014). 
Se denomina K a la matriz de rigidez en coordenadas principales; para el inicio del 
calculo cuando trabaja la estructura en el rango elastico, esta matriz resulta. 
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KL= 



Para estructura con disipadores ADAS: 

17831.33 -2197.03 15.79 -83.76 14.85 14.73 -1.38 

-2197.03 28447.70 -14930.6 -7477.99 690.11 2225.43 -309.16 

15.79 -14930.6 15637.41 -24.60 -663.41 -7.20 -41.82 

-83.76 -7477.99 -24.60 24066.60 -14262.4 -6031.31 1604.47 

14.85 690.11 -663.41 -14262.4 14972.19 -658.52 -107.33 

14.73 2225.43 -7.20 -6031.31 -658.52 8114.36 -3171.48 

-1.38 -309.16 -41.82 1604.47 -107.33 -3171.48 1949.59 

Para estructura con disipadores TADAS: 

18895.17 -3276.92 23.62 -125.22 22.10 22.12 -2.06 

-3276.92 29544.39 -14929.1 -7445.70 677.55 2222.90 -308.30 

23.62 -14929.1 15979.36 -369.39 -656.22 -20.56 -39.80 

-125.22 -7445.70 -369.39 24415.80 -14260.5 -6027.47 1602.26 

22.10 677.55 -656.22 -14260.5 15319.28 -1002.87 -114.59 

22.12 2222.90 -20.56 -6027.47 -1002.87 8456.41 -3164.32 

-2.06 -308.30 -39.80 1602.26 -114.59 -3164.32 1949.75 



3.2 Matriz de masas 



KL= 



Sea m 1 la masa del primer piso; m 2 ,m 3 ,m 4 , son las masas de los pisos 2 a 
4, las mismas que se calculan en base a las cargas indicadas en la figura 15. En 
este caso, todas las masas son iguales. La matriz de masas se obtiene evaluando 
la energia cinetica de la estructura, T c 

T c = |mu 2 + - J 9 2 (19) 

Donde m es la masa; v es la velocidad lineal evaluada en el centro de 
masas; J es el momento de inercia de la masa evaluada en el centro de masas; 8 
es la velocidad angular. (Aguiar 2012) En el modelo se desprecia J con lo que la 
matriz de masas en T*s 2 /m resulta. 
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3.3 Propiedades dinamicas 

Notese que la matriz de masas tiene ceros en la diagonal, por lo que 
interesa ver los periodos de vibracion y los respectivos modos de vibracion, los 
mismos que se indican en la tabla 1, destacando que los modos estan 
normalizados de tal manera que en el tope de la estructura sea igual a 1 . 

Tabla 1 Periodos y modos de vibracion. 



Para estructura con disipadores ADAS: 



MODOS 
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Para estructura con disipadores TADAS: 


MODOS 


1 


2 


3 


4 


PERIODOS 


0.36798 


0.12465 


0.06024 


0.03898 
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Los modos de vibracion colocados son los pertenecientes a cuando la 
masa modal efectiva es igual al 100%, y este numero de modos corresponded a 
los pisos donde existe masa diferente de cero. 

Finalmente, interesa conocer los factores de participacion modal r m 
debido a que se necesitan para explicar el pushover multimodal. 

= 0® t Mb ( 20 ) 

0(0* m 0(0 

Donde 0® es el modo de vibracion i; M es la matriz de masas; b es el 
vector de incidencia de las coordenadas principales con respecto al movimiento 
del suelo. Cuando se trabaja con coordenadas absolutas, como en este caso el 
vector b es unitario y para el presente caso, esta compuesta por 7 unos. 
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Los factores de participacion indican que tanto participa el modo en la 
respuesta. En la tabla 2 se indican los factores de participacion y se observa que 
el factor de participacion del primer modo es el mas importante y esto se refleja en 
la curva de capacidad sismica resistente, como se vera mas adelante. 



Tabla 2 Factores de participacion modal. 



Modo 


Factor Participacion 
ADAS 


Factor de Participacion 
TADAS 


1 


1 .32371 


1 .33291 


2 


-0.43725 


-0.44158 


3 


0.15021 


0.14386 


4 


-0.03667 


-0.03519 



3.4 Programas de CEINCI-LAB que encuentran matriz de rigidez, 
masas y propiedades dinamicas. 



A continuacion se presentan los programas del sistema de computacion 
CEINCI-LAB que obtienen exclusivamente la matriz de rigidez, matriz de masas y 
las propiedades dinamicas, que se ha presentado en este apartado. 



Tabla 3 Descripcion de programas 



Programa 


Descripcion 


Geom 


Determina las propiedades de las 
secciones, como el area, el momento de 
inercia, el momento de fluencia, la rigidez 
axial de fluencia y la rigidez de fluencia. 


Structure_Geom_2D_v2 


Programa que determina la geometna de 
la estructura, los grados de libertad, la 
matriz que contiene los vectores de 
colocacion, la longitud, seno y coseno de 
los elementos. 


espe_cap3 


Programa que determina las propiedades 
dinamicas de la estructura como periodos, 
modos de vibracion, factores de 
participacion modal, masa modal efectiva, 
condensacion de matriz de rigidez, etc. 


lteration_Vy 


Sirve para determinar en forma exacta el 
punto de fluencia en una seccion, cuando 
se pasa de la zona elastica a la plastica, 
regresa y disminuye el incremento de 
carga. 


Dibujo_Din2 


Con este programa se tiene la opcion de 
indicar cada cuantos ciclos de carga dibuja 
la estructura y coloque con colores en los 
elementos la zona del diagrama momento 
curvatura en que se halla. El color verde es 
para la zona post fluencia y el color rojo 
para la zona residual. 


masas2 


Encuentra la matriz de masas en el modelo 
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en que se considera como un elemento al 
disipador. Las masas actuan en el piso 
pero al considerar un elemento al 
disipador, en la parte inferior no hay masa 
por lo que su valor es cero. 


fuerzas_0 


Encuentra las fuerzas en coordenadas 
locales, en cada ciclo de carga. Luego son 
almacenadas en la variable Tabla2. 


Momento_Rotacion_ASCE41 B 


Si momento es positivo en vigas, As sup se 
encuentra en la parte superior y As inf en la 
parte inferior; caso contrario se invierte la 
armadura. Ademas llama a programas que 
determinan los parametros de ASCE 41 


Momento_Rotacion_ASCE41 C 


Similar al programa anterior pero para 
columnas, donde interviene la carga axial. 


Momento_Rotacion_ConcrASCE41 B 


Determina los parametros de ASCE 41 
para el caso de vigas de hormigon armado. 


Momento_Rotacion_ConcrASCE41 C 


Determina los parametros de ASCE 41 
para el caso de columnas de hormigon 
armado. 


Momento_Rotacion_AceroASCE41 B 


Determina los parametros de ASCE 41 
para el caso de vigas de acero. 


Momento_Rotacion_AceroASCE41 C 


Determina los parametros de ASCE 41 
para el caso de columnas de acero. 


krigdez_ASCE41-2Mfi 


Pasa de la relacion momento-rotacion a 
momento-curvatura y encuentra rigideces a 
flexion en nudo inicial, centro de luz y nudo 
final con las cuales halla la matriz de 
rigidez de elementos vigas y columnas 
usando el modelo incluido en K1_eqcorte4 
(Mora y Aguiar (2015)) y el modelo de 
Giberson, luego ensambla la matriz de 
rigidez total de la estructura. 


krigidez_braces_eq 


Determina la rigidez axial equivalente de 
contraviento-disipador. Ademas puede 
calcular la rigidez axial solo de diagonales 
articuladas. 


krigidez_damper2 


Obtiene la rigidez de disipador con modelo 
2 A y ensambla en la matriz de rigidez de 
la estructura. 


krigidez_damper3 


Similar al programa anterior pero con 
modelo 2 B. 


PlasticHinge 


Determina la longitud plastica. 


AxialTADAS 


Para el modelo bilineal fuerza- 
desplazamiento del elemento disipador 
TADAS identifica la rama en que se 
encuentra para la rigidez axial. 


AxialADAS 


Similar al anterior pero con disipador 
ADAS. 
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3.5 Programa desarrollado que obtiene matrices de rigidez, masa y 
propiedades dinamicas 

• Programa en Matlab para hallar la Matriz de Rigidez de una 
estructura sin arriostrar: 

krigdez_ASCE41_2Mfi 

function 

[SS, icon, iconu, iconc,Table3,K3, Hinge, tetay, TableMf i] =krigidez 
_ASCE41_2Mf i (ngl, Tab lei, Table2, Table3, Table4, L, seno, coseno, VC 
, icon, iconu, iconc, j j , Pa, pq, Q2cg, Hinge, tetay) 

mbr=length (L) ; SS=zeros (ngl) ; TableMf i=zeros (mbr, 12) ; 

K3=zeros (6, 6, mbr) ; 

STl=size (Tablel, 2) /3; 

for i=l :mbr 

if Tablel (i, 1) ~=1 . 1 && Tablel ( i , 1 ) ~=1 . 2 

tetayi=tetay(i, 1) ;tetayj=tetay(i,2) ; 

lon=L ( i ) ; sen=seno ( i ) ; cosen=coseno ( i ) ; 

tetai=abs (Table4 (i, 1) ) ; tetaj=abs ( (Table4 (i, 2) ) ) ; 

Mi=-Table2 (i, 3) ;Mj=Table2 (i, 6) ; 

Vi=abs (Table2 (i, 2) ) ;Vj=abs (Table2 (i, 5) ) ; 
if (sen>=- . 707106 && sen<=. 707106) 

[ai, bi, ci, tetayi, EAa, Elya, Myi, pa] =Momento_Rotacion_ASCE4 IB (Ta 
blel (i, 1 : ST1) , Ion, Mi, tetai, tetayi, Vi) ; 

[~,~,~, ~ i EAo, Elyo, ~] =Momento_Rotacion_ASCE41B (Tablel (i, ST1 + 1 : 
ST1*2) , Ion, 0,0,0,0); 

[aj , bj , cj , tetayj , EAb, Elyb, Myj , pb] =Momento_Rotacion_ASCE4 IB (Ta 
blel (i,STl*2+l:STl*3) , Ion, Mj , tetaj , tetayj , Vj ) ; 

else 

[ai, bi, ci, tetayi, EAa, Elya, Myi, pa] =Momento_Rotacion_ASCE4 1C (Ta 
blel (i, 1 : ST1) , Ion, Mi, Pa (i) , tetai , tetayi , Vi ) ; 
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[~,~,~, ~, EAo, Elyo, ~] =Momento_Rotacion_ASCE41C (Tablel (i, ST1 + 1 : 
ST1*2) , Ion, 0, Pa (i), 0,0,0); 

[aj , bj , cj , tetayj , EAb, Elyb, Myj , pb] =Momento_Rotacion_ASCE4 1C (Ta 
blel (i,STl*2 + l:STl*3) , Ion, Mj , Pa (i) , tetaj , tetayj , Vj ) ; 

end 

tetay (i, : ) = [tetayi, tetayj ] ; 

g, 
o 

%i linear y j linear 

EIa=EIya; EIo=EIyo; EIb=EIyb; 

g, 
o 

Mreal= [Mi,Mj ] ; 

[Lpa, Lpb, beta_a, beta_o, beta_b] =PlasticHinge (Tablel (i, :) , Ion, M 
real, Ela, EIo, Elb, Myi, Myj ) ; 

Lo=lon-Lpa-Lpb; 

Mui= (ai/Lpa) *EIya*pa+Myi ; 

Muj= (aj /Lpb) *EIyb*pb+Myj ; 

%i nonlinear 

if abs (Mi) >=Myi 

EIa=pa*EIya; 

if Table3 (i, 1) ==0 

icon=icon+l ; 

Table3 (i, 7) =icon; Table3 (i, 1) =j j ; 

end 
end 

%j nonlinear 

if abs (Mj ) >=Myj 

EIb=pb*EIyb; 

if Table3 (i, 2) ==0 

icon=icon+l ; 

Table3 (i, 8 ) =icon; Table3 (i,2)=j j; 
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end 
end 

O 

if (Table3 (i, 5) ~=0 || Table3 ( i , 6 ) ~=0 ) || (tetai>tetayi+bi | 
teta j >tetayj +b j ) 

if tetai>tetayi+bi || Table3 ( i , 5 ) ~=0 

if Table3 (i, 5) ==0, 
iconu=iconu+l ; Table3 (i, 11) =iconu; Table3(i,5)=jj; 
end 

T12=[0 1/lon 0 0 -1/lon 1; 
-1 0 0 1 0 0]; 

Kl=[ (6*EIb*EIo) / (ion* (Elb+EIo) ) 0 
0 ( E Aa * E Ab * E Ao ) / (EAa*EAb*Lo+EAa*EAo*Lpb+EAb*EAo*Lpa) ] ; 

end 

if tetaj>tetayj+bj || Table3 ( i , 6 ) ~=0 

if Table3 (i, 6) ==0, 
iconu=iconu+l ; Table3 (i, 12) =iconu; Table3(i,6)=jj; 
end 

T12=[0 1/lon 1 0 -1/lon 0 ; 

-1 0 0 1 0 0]; 

Kl=[ (6*EIa*EIo) / (Ion* (Ela+EIo) ) 0 
0 ( E Aa * E Ab * E Ao ) / ( EAa * EAb * Lo+EAa * EAo * Lpb+EAb * EAo * Lpa ) ] ; 

end 

if tetai>tetayi+bi && teta j >tetayj +b j 

Kl= (EAa*EAb*EAo) / ( EAa * EAb * Lo+EAa * EAo * Lpb+EAb * EAo * Lpa ) ; 
T12= [-1 0 0 1 0 0]; 
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end 

else 

T12=[0 1/lon 1 0 -1/lon 0; 
0 1/lon 0 0 -1/lon 1; 

-1 0 0 1 0 0]; 

[Kl] =Kl_eqcorte4 (EIya,EIa,EIo,EIyb,EIb, EAa , EAa , EAo , EAb , EAb , ma 
x (lon/3-Lpa, 0) , Lpa, lon/3+min (lon/3-Lpa, 0) +min (0, Ion/ 3- 
Lpb) , Lpb, max ( 0 , lon/3-Lpb) , beta_a) ; 

end 

if (Mi>Mui && tetai<=tetayi+bi) | | (Mj>Muj && 
tetaj<=tetayj+bj ) || Table3 ( i , 3 ) ~=0 || Table3 ( i , 4 ) ~=0 

%i rotula pi 

if Mi>Mui && tetai<tetayi+bi || (Table3 (i, 3) ~=0 && 
Table3 (i, 5) ==0) 

Kl (1, 1) =- (Q2cg (i, 3) *lon-ci*Myi*lon) / (cosen*pq (2, i) - 
cosen*pq (5, i) +lon*pq (3, i) -pq (1, i) *sen+pq (4, i) *sen) ; 

Kl (1,2)=0;K1 (2, 1)=0; 

if length (Kl) ==3 

Kl (2, 2) = (6*EIb*EIo) / (Ion* (Elb+EIo) ) * (1/2* (2+beta_o) / (l+2*beta 
_o) ) ; 

end 

if Table3 (i, 3) ==0 
iconc=iconc+l ; 

Table3 (i, 9) =iconc+l; Table3 (i, 3) =jj ; 

end 
end 

%j rotula pi 

if Mj>Muj && tetaj<tetayj+bj || (Table3 (i, 4) ~=0 && 
Table3 (i, 6) ==0) || (Table3 (i, 4) ~=0 && Table3 ( i , 6 ) ==0 ) 

KlcMy=- (Q2cg (i, 6) *lon-cj *Myj *lon) / (cosen*pq (2, i) - 
cosen*pq ( 5 , i ) +lon*pq ( 6 , i ) -pq ( 1 , i ) *sen+pq ( 4 , i ) *sen) ; 
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Kl (1,2)=0;K1 (2, 1)=0; 
if length (Kl) ==3 

Kl (1, l) = (6*EIa*EIo) / (Ion* (Ela+EIo) ) * (1/2* (2+beta_o) / (l+2*beta 
_o) ) ; 

Kl (2, 2) =KlcMy; 

else 

Kl (1, 1) =KlcMy; 

end 

if Table3 (i, 4) ==0 
iconc=iconc+l ; 

Table3 (i, 10) =iconc+l; Table3 (i, 4) =j j ; 

end 
end 

end 



K2=T12 ' *K1*T12; 

T23= [cosen sen 0 0 0 0 

-sen cosen 0 0 0 0 

0 0 1 0 0 0 

0 0 0 cosen sen 0 

0 0 0 -sen cosen 0 

0 0 0 0 0 1 ] ; 

K3 ( : , : , i) =T23 ' *K2*T23; 

for j =1 : 6 

J=VC(i, j) ; 

if J==0 

continue 

end 

for m=l : 6 

mm= VC ( i , m ) ; 



if mm==0 
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continue 

end 

SS (J, mm) =SS (J, mm) +K3 ( j ,m, i) ; 

end 

end 

f iyi=Myi/EIya; f iyj=Myj /Elyb; f iui= (Mui- 

Myi) /pa/EIya+f iyi; f iuj= (Muj -Myj ) /pb/EIyb+f iyj ; 

if jj==l 

f iui=min (10*f iyi, (Mui- 

Myi) /pa/EIya+f iyi) ; fiuj=min (10*fiyi, (Muj-Myj ) /pb/EIyb+f iyj ) ; 
Mui=min(Mui, (f iui-f iyi) *pa*EIya+Myi) ;Muj=min (Muj , (fiuj- 
f iyj ) *pb*EIyb+Myj ) ; 

End 

TableMfi (i, : ) = [0, f iyi, f iui, 0,Myi,Mui, 0, fiyj , fiuj , 0,Myj , Muj ] ; 

end 

end 

• Programs en Matlab para obtener la Matriz de Masas: 
masas2 

function [Mt, M] =masas2 (Fm, NI , NJ, Y, L, np) 
M=zeros (np) ; 

Niv=unique ( [Y (NI) , Y (NJ) ] ) ; 
for oo=l : np 

for mm=l : size (Fm, 1 ) 

if Y (NI (Fm (mm, 1) ) ) ==Niv (oo+l) 

M (oo, oo) =M (oo, oo) +L (Fm (mm, 1) ) *abs (Fm (mm, 2) ) / 9 . 8 ; 

end 
end 



end 
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Mt=sum (sum (M) ) ; 

• Programa en Matlab para obtener propiedades dinamicas: 
espe_cap3 

function 

[FI, T, gama, alfa, OM, KL] =espe cap3 (K,M,Mt, np, ngl, nmodos) 

kaa=K (1 :np, 1 :np) ; kab=K (1 :np, np+1 :ngl) ; kbb=K (np+1 :ngl, np+1 :ngl 
) ; kba=kab ' ; 

KL=kaa-kab*kbb A -l*kba; 
[T,FI,OM]=orden_eig(KL,M) ; 
gama=zeros (nmodos, 1) ; alfa=gama; 

FI=FI . /repmat (FI (np, : ) , np, 1 ) ; % Vibration Modes normalized to 
1 on the 1st floor 

J=ones (np, 1 ) ; 

for i=l :min (size (FI, 2) , nmodos) 
fi_n=FI ( : , i) ; 

Ln= (f i_n ' *M*J) ;Mn= (f i_n ' *M*f i_n) ; 

gama ( i ) =Ln/Mn; %gama=Ln/Mn Factor de participacion modal 

alfa (i) = (f i_n ' *M*J) A 2/ (Mt* (fi_n ' *M*f i_n) ) ; %alf a=Mn_a/Mt; 
sn=gama*M*fi n 

end 



4 PUSHOVER MULTIMODAL 

En los primeros trabajos del Pushover se remontan a la decada de los 
anos setenta y ochenta del siglo pasado, las cargas laterales se aplicaban en 
forma monotonica en cada piso de la estructura, las mismas que eran obtenidas 
en base al primer modo de vibracion (Freeman, 1975; Freman et at. 1978). Por 
esos tiempos el desarrollo informatico era incipiente por lo que una buena parte de 
los Proyectistas estructurales o los investigadores no le dieron mayor importancia 
al asunto. 
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Con el gran desarrollo informatico que se tiene en la ultima decada del 
siglo XX, los investigadores empiezan a estudiar lo publicado y empiezan a 
cuestionar el Pushover se dan cuenta que en estructuras irregulares el aplicar 
cargas laterales en funcion del primer modo de vibracion no es adecuado y 
recomiendan que se o haga en base a todos los modos de vibracion. Chopra y 
Goel (1999, 2001, 2002) 

Es asi como nacio el Pushover multimodal que se presenta en este 
articulo, destacando que actualmente se esta trabajando con el Analisis 
Incremental Dinamico IDA, en el cual las fuerzas laterales que se aplican a la 
estructura pueden ser proporcionales a la forma espectral del primer modo, como 
lo recomienda el ASCE 7-05. (Vielma y Cando, 2014; Vielma 2015). Con esta 
pequeha introduccion se indica a continuacion el Pushover multimodal. 

En un analisis modal se incrementa el cortante basal en cada ciclo de 
carga y este se lo reparte para cada modo de vibracion usando una distribucion de 
fuerzas laterales s* = m4> h como se indica en Chopra (2014), siendo m la masa 
de cada piso, <fo es el modo de vibracion i. 

Luego la reparticion de la carga lateral por piso se la hace de acuerdo a la 
participacion modal r t ya indicada en el apartado anterior, donde i representa el 
modo de vibracion. 

Entonces, se obtiene un vector de cargas laterales para cada ciclo y para cada 
modo de vibracion igual a: 

Qa n = npis T 9iM>K T n *AV (21) 

Donde: diag(M) es el vector de la diagonal de la matriz de masas de la 
estructura (si la matriz es diagonal); (p t es el modo de vibracion i; T t es el factor de 
participacion modal; AV es el incremento del cortante basal. 

Es conveniente que el valor AV sea pequeho pero va a demandar mucho 
de tiempo de ejecucion y se debe tener presente que los programas realizan un 
ajuste cuando se pasa del rango elastico al inelastico, ajuste que consiste en 
regresar al ciclo de carga anterior y reducir el incremento de carga hasta llegar al 
punto en que se pasa al rango no lineal. 

4.1 Criterio de combinacion modal CQC 



Si se realiza un analisis paso a paso de historia en el tiempo y separamos 
los resultados obtenidos para cada modo de vibracion de la estructura se podra 
observar que en el transcurso de la respuesta de las estructuras ante un 
acelerograma habra picos absolutos maximos en diferentes sitios, esto hace algo 
dificil combinar exactamente los resultados ya que un analisis pushover representa 
el pico absoluto maximo de respuesta de un analisis dinamico en historia en el 
tiempo. 
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En Aguiar (2012) se presenta unos diez criterios de combinacion modal lo 
que demuestra que el tema es bastante complejo. De estos criterios, el mas 
utilizado y acogido en algunas normativas sismicas es el CQC (Complete 
Quadratic Combination, Chopra, 2001), en que toma en cuenta la posibilidad de 
acoplamiento entre los modos. 

El CQC realiza una combinacion de suma de cuadrados mas un termino 
que combina las respuestas con un factor que dependera de la frecuencia y factor 
de amortiguamiento de cada modo de vibracion, este termino puede ser positivo o 
negative Se presenta a continuacion las formulas utilizadas para combinar las 
respuestas r con el criterio CQC: 

N N 



r " =Z2>; 1 r } (22) 



=i ./=i 



Lo que es igual a escribirla de la siguiente forma, la cual resume lo descrito 
en el parrafo anterior: 



r = (xU r n 2 + ^i^U Pim) 



!/2 ( 23 ) 



i*j 



Donde, 



P 'J = 1 h v a = ^T (24)y(25) 

(l-a 2 ) + 4£ 2 a(l + a) 2 W ni 

P = ( 26 ) 

' J (l-a 2 ) 2 + 4^.a(l + « 2 )+4fe 2 + ^ 2 y 

Donde iyj corresponden al i-esimo y j - esimo modo de vibracion; 
W nU es la frecuencia de vibracion en el modo i; W nj es la frecuencia de vibracion 
en el modo son los factores de amortiguamiento en los modos de vibracion i 

yj- 



4.2 Factores de amortiguamiento 

Para la combinacion modal vista en el sub apartado anterior, se necesita el 
factor de amortiguamiento de la estructura en cada modo de vibracion. Los 
factores utilizados para los modos diferentes del primer modo de vibracion se los 
tomo de Naoki et al. (2003) que se encuentran detallados en Aguiar (2012). 

Se tomo el criterio de factores de amortiguamiento de Naoki et al. (2003), 
debido a su sustento practico, y este se resume en las siguientes formulas: 
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• Para estructuras de Hormigon Armado. 

6 = 1-4 * ( 27 ) 

• Para estructuras de acero. 

f, = (1.3 - 1.4) * ( 28 ) 

• Para estructuras de Mixtas 

& = (1.7 - 1.8) * ( 29 ) 

Donde £ es el factor de amortiguamiento en el modo i; es el modo de 
vibracion en el modo i - 1 . Para el primer modo de vibracion se trabajo con el 
amortiguamiento intnnseco de la estructura igual a 0.05 y a partir de este valor se 
obtuvieron los restantes factores de amortiguamiento. 



5 ESTRUCTURA DE HORMIGON ARMADO 



Dos estructuras se han analizado, una de hormigon armado de 4 pisos y 
otra de acero de 6 pisos. En cada una de ellas se ha contemplado la colocacion de 
dos tipos de disipadores a saber: ADAS y TADAS. De tal manera que se tienen 4 
casos estudiados. 

La estructura de hormigon armado, es la indicada en las figuras 13 a 15 y 
para explicar un poco mas la forma como se obtienen las cargas incrementales 
que se aplican en cada modo de vibracion, se indica en las tablas 4 y 5 las cargas 
aplicadas en (T), para el rango elastico para la estructuras con ADAS y TADAS, 
respectivamente. 

Tabla 4 Cargas laterales en (T), aplicadas en el rango elastico, para cada ciclo de 
carga. Estructura con disipadores ADAS 





1 


2 


3 


4 


1 


0.0000000 


0.0000000 


0.0000000 


0.0000000 


2 


0.0021037 


-0.0046398 


0.0045037 


-0.0017835 


3 


0.0000000 


0.0000000 


0.0000000 


0.0000000 


4 


0.0062178 


-0.0088846 


0.0013665 


0.0017452 


5 


0.0000000 


0.0000000 


0.0000000 


0.0000000 


6 


0.0101430 


-0.0040760 


-0.0043755 


-0.0011412 


7 


0.0131943 


0.0075375 


0.0019570 


0.0003316 
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Tabla 5 Cargas laterales en (T), aplicadas en el rango elastico para cada ciclo de 
carga. Estructura con disipadores TADAS 





1 


2 


3 


4 


1 


0.0000000 


0.0000000 


0.0000000 


0.0000000 


2 


0.0018888 


-0.0041283 


0.0039616 


-0.0015908 


3 


0.0000000 


0.0000000 


0.0000000 


0.0000000 


4 


0.0057064 


-0.0082885 


0.0013953 


0.0015229 


5 


0.0000000 


0.0000000 


0.0000000 


0.0000000 


6 


0.0094615 


-0.0042100 


-0.0038509 


-0.0009709 


7 


0.0127411 


0.0070556 


0.0016158 


0.0002704 



En la figura 16 se presenta la curva de capacidad sismica resistente, 
hallada en cada modo de vibracion, de la estructura con disipadores ADAS. En la 
parte superior derecha se encuentra la contribucion del primer modo, en esta 
grafica la escala vertical Mega a 35 T; a la derecha se tiene la contribucion del 
segundo modo, aqui la escala Mega a 9 T, practicamente 4 veces menor; en la 
parte inferior se indica la contribucion del tercer y cuarto modo sus escalas son 
1 .63 T., y 0.7 T. Esto indica que en la presente estructura es el primer modo el que 
mas aporta a la respuesta estructural. 

Para encontrar la curva de capacidad sismica resistente total se aplica el 
criterio de combinacion CQC, encontrando la grafica que se indica en la figura 17, 
con linea continua. 

En la figura 17, con linea entrecortada se observa, con bastante dificultad 
la contribucion de cada uno de los modos; lo que se aprecia es la contribucion del 
primer modo, de los restantes modos se halla en la parte inferior. 

Se obtuvo la curva de capacidad sismica, mediante un pushover 
monotonico, en el cual las cargas laterales se aplican proporcionalmente de 
acuerdo al primer modo de vibracion y se encontro que la respuesta difiere de la 
curva que se obtiene con el pushover multimodal, como se observa en la figura 18, 
donde se encuentran los resultados de la figura 17 mas la nueva curva de 
capacidad que se halla con el pushover monotonico. 
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Desplazamiento en tope [m] 

3er Modo de Vibracion 



Desplazamiento en tope [m] 

4to Modo de Vibracion 



Figura 16 Curvas de Capacidad con modos de vibracion de Estructura de 
Hormigon de 4 pisos con disipadores ADAS, con elemento disipador. 




0.1 0.2 0 3 0 4 

Desplazamiento en tope [m] 



0 5 



Figura 17 Curvas de Capacidad resultante de Estructura de Hormigon de 4 pisos 
con disipadores ADAS, con elemento disipador 
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Figura 18 Curva de capacidad sismica encontrada en estructura de hormigon 
armado, de cuatro pisos con disipadores ADAS. Con dos modelos de calculo del 
Pushover: Monotonico y Multimodal. 

Antes de comparar las curvas de capacidad encontradas con el Pushover 
multimodal y el Pushover monotonico, es importante destacar que el cortante 
basal, que se tiene en el eje vertical de las curvas de capacidad se obtiene 
sumando el incremento de cortante, AV que se va aplicando a la estructura y que 
se indico en el apartado 4. El criterio de combinacion modal CQC, se aplica en los 
desplazamientos, momentos y fuerzas que se obtienen en cada ciclo de carga. 

En efecto para un determinado ciclo de carga, se tiene tantos estados de 
carga lateral que se aplican a la estructura, como modos de vibracion se tengan; 
de tal manera que para cada iteracion se tiene una respuesta en desplazamientos, 
fuerzas y momentos, en cada modo de vibracion y es ahi cuando se aplica el CQC 
y se hallan la resultante. 

Al cambiar los disipadores ADAS por los TADAS, se obtuvo resultados 
similares por lo que se omite su presentacion y unicamente se presenta la figura 
19, que corresponde a la curva de capacidad sismica que se obtiene con el 
pushover monotonico y con el pushover multimodal. En esta grafica tambien se 
incluye con lineas entrecortadas la contribucion de cada uno de los modos; los 
comentarios son similares a los realizados en el parrafo anterior. 
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Desplazamiento en tope [m] 

Figura 19 Curva de capacidad sismica encontrada en estructura de hormigon 
armado, de cuatro pisos, con disipadores TADAS. Empleando dos modelos de 
analisis: Monotonico y Multimodal. 



6 ESTRUCTURA DE ACERO 

En la figura 20 se presenta la estructura de acero de 6 pisos y tres vanos, 
con disipadores de energia colocados en el vano central, en todos los pisos. En 
Aguiar et al. (2015, 2) se indica la forma como se halla el diagrama momento 
rotacion de estos perfiles utilizando las recomendaciones del ASCE 41 de 2011 y 
posteriormente como se pasa al diagrama momento curvatura en funcion de la 
longitud plastica. Por cierto, en los programas que obtienen el diagrama momento 
curvatura, indicados en la tabla 3, para el rango elastico solo encuentran el 
momento de fluencia y curvatura de fluencia; cuando ingresan al rango no lineal 
ahi empiezan a calcular el punto U, con los coeficientes del ASCE 41 . 

En la figura 21 se indica a la izquierda los disipadores de energia tipo 
TADAS y a la derecha los ADAS. De tal manera que se analizaron dos casos con 
cada uno de estos disipadores. 

En la figura 22 se presentan las curvas de capacidad sismica, encontradas 
en cada uno de los modos de vibracion, en la estructura con disipadores ADAS. 
En la parte superior se hallan las curvas para los tres primeros modos de vibracion 
y en la parte inferior para tres ultimos modos. Para el primer modo, la escala 
vertical que indica el cortante basal Mega hasta las 400 T (curva de capacidad) en 
cambio para el segundo hasta las 120 T y en el tercer modo el cortante es del 
orden de las 70 T. De tal manera que el primer modo es fundamental en la 
respuesta de la estructura, en este caso en la curva de capacidad sismica 
resistente. 
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Figura 20 Estructura de acero con disipadores de energia colocados sobre 
diagonales en forma de "V" invertida. 
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Figura 21 Geometna de los disipadores de energia ADAS, a la izquierda y TADAS 

a la derecha. 
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Figura 22 Curvas de Capacidad con modos de vibracion de Estructura de Acero de 
6 pisos con disipadores ADAS, con elemento disipador. 

La escala vertical, para el cortante basal del cuarto modo Mega a 20 T., del 
quinto a 1 .4 T y del sexto a 0.06 T. En las figuras se esta indicando con elemento 
disipador, esto significa que el disipador de energia fue modelado como un 
elemento mas de la estructura, en contraste con el modelo de la diagonal 
equivalente en el cual la diagonal de acero y el disipador se modelan como una 
sola diagonal equivalente. 

En la figura 23 se presenta con Ifnea entrecortada, la contribucion de cada 
uno de los modos y con linea continua la resultante que se obtiene al aplicar el 
criterio de combinacion modal CQC. En este caso se observa que la resultante se 
aleja de la curva del primer modo y que tiene valores muy inferiores para 
desplazamientos altos. 

En la figura 24 se compara las curvas de capacidad sismica que se hallan 
con dos modelos numericos de calculo del Pushover, el uno el monotonico que 
trabaja solo con el primer modo de vibracion y el segundo el multimodal 
considerando todos los modos de vibracion. En esta grafica se nota una mejor 
coincidencia de la curva de capacidad hallada con el modelo monotonico, con la 
curva que se encuentra en el primer modo de vibracion. 

Otro aspecto interesante de la figura 24 es que las curvas que se 
encuentran con los dos modelos: monotonico y multimodal, son bastante 
diferentes. 
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Figura 23 Curvas de capacidad sismica resultante de Estructura de Acero de 6 
pisos con disipadores ADAS, con elemento disipador. Con linea entrecortada se 
indica la contribucion de cada modo y con continua la resultante. 
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Figura 24 Curvas de Capacidad Modal y Monotonica de Estructura de Acero de 6 
pisos con disipadores ADAS, con elemento disipador y contribucion de cada modo 

de vibracion. 
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Nuevamente cuando se cambia de disipador de ADAS a TADAS, el 
comportamiento que se observa es muy similar, razon por la cual solo se presenta 
la curva de capacidad sismica que se halla con los dos modelos: monotonico y 
multimodal, en la figura 25. Se observa una muy buena aproximacion en los 
resultados hallados con estos dos modelos. 




Figura 25 Curvas de capacidad sismica resistente, con dos modelos numericos de 
calculo: Modal y Monotonica, de una estructura de acero de 6 pisos con 
disipadores TADAS, considerando elemento disipador. 



7 COMENTARIOS Y CONCLUSIONES 

Se ha presentado con bastante detalle, la forma de calculo de la curva de 
capacidad sismica resistente de una estructura empleando la tecnica denominada 
Pushover Multimodal. A sabiendas de que existen muchas publicaciones al 
respecto, pero una buena parte de ellas se limitan a utilizar determinados 
programas de ordenador y en otras la explicacion es bastante general. Se aspira 
con esta publicacion haber ilustrado el calculo del Pushover empleando todos los 
modos de vibracion para aplicar las cargas laterales y aplicando el criterio de la 
combinacion cuadratica completa CQC en los desplazamientos, fuerzas y 
momentos que se generan en cada ciclo de carga. 

Por otra parte se ha encontrado la curva de capacidad sismica resistente 
en una estructura con disipadores de energia ADAS o TADAS, los mismos que se 
hallan sobre contravientos Chevron. Los disipadores han sido modelados como un 
elemento mas de la estructura y la matriz de rigidez fue obtenida considerando 
cinco dovelas en el disipador. En los restantes elementos tambien se obtuvo la 
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matriz de rigidez en coordenadas globales considerando tres dovelas. De tal 
manera que en las columnas y vigas se esta trabajando con un modelo de 
plasticidad que toma en cuenta la longitud del elemento que ingresa al rango no 
lineal. 

Los resultados encontrados del analisis Pushover Multimodal han sido 
comparados con los que se obtienen con el Pushover Monotonico, encontrando 
una buena aproximacion, debido a que la estructura es simetrica en elevacion. 

Finalmente todos los programas que se han desarrollado han sido 
incorporados al sistema de computacion CEINCI-LAB. 
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